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ABSTRACT 


Two  models  of  atmospheric  flow,  "compensation"  and  "moist 
over-compensation"  are  found  to  be  useful  in  the  study  of  tropical  weather 
in  conventional  and  satellite  observations.  A  numerical  weather  predic¬ 
tion  model  simulating  important  physical  processes  in  the  low  latitudes  is 
developed  to  test  the  feasibility  of  using  ATS  information  to  infer  future 
weather  distribution.  Diagnostic  studies  using  ATS  films  and  conven¬ 
tional  synoptic  weather  maps  describe  the  weather  situations  in  the  Cen¬ 
tral  North  Pacific  in  terms  of  the  theoretical  models.  Results  of  the  pre¬ 
diction  experiment  and  future  use  of  ATS  type  data  in  numerical  prediction 
models  are  discussed  and  application  recommended.  Appendix  A  gives  a 
description  of  one  method  of  following  cloud  motion  in  satellite  photographs. 
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I.  FORECAST  METHODS  FROM  SATELLITE  PHOTOGRAPHS 


John  C.  Freeman  and  Sam  Rosenberg 


A.  Introduction 

The  objective  of  this  study  ls  to  develop  a  numerical  weather 
Di-edictlon  model  simulating  important  physical  processes  in  the  tropical 
and  subtropical  atmosphere  which  can  be  applied  to  data  furnished  by 
ATS  satellite  photographs,  in  order  to  test  the  feasibility  of  using  ATS 
Information  to  infer  future  weather  distribution. 

Parallel  with  the  theoretical  development  of  the  prediction  model, 
it  was  necessary  to  carry  out  preliminary  diagnostic  studies,  using  ATS 
films  and  microfilm  copies  of  manuscript  maps  from  12  March  to  31  May 
1967.  Only  hemispheric  NMC  surface  charts  were  available  for  the  month 
of  March,  1967.  For  the  remainder  of  the  period  studied  the  maps  prepared 
by  the  ESSA  Weather  Bureau  in  Honolulu  were  used.  Analyses  were  made 
available  for  the  surface,  700,  500  and  250  millibar  levels. 

The  perusal  of  t.nese  data  led  to  the  selection  of  a  sequence  of 
weather  (Section  II)  for  the  application  of  the  numerical  model  described 
in  Section  I.  The  April  case  selected  for  the  series  of  prediction  experi¬ 
ments  in  Section  III  was  one  which  is  frequently  observed  over  the 
tropical  and  subtropical  oceans,  with  a  typical  broad  zone  of  tropical 
easterlies.  The  series  of  ATS  photos  showed  an  increase  in  cloudiness  and 
change  in  circulation  during  the  selected  forecast  period.  The  March  case 
study  presented  in  Section  IV  was  selected  for  analysis  because  it  was 
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representative  of  several  good  examples  where  the  information  contained 
in  the  sequences  of  ATS  photos  significantly  augmented  conventional  data. 

In  addition,  examination  of  these  data  has  provided  some  insight 
for  possible  future  application  of  this  model  into  ATS  type  data  for  the 
numerical  prediction  of  clouds,  precipitation  and  circulation.  (Section  IV) 
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R,  Compensation  and  Clouds 

The  concept  of  compensation  is  offered  as  a  key  to  understanding 
and  predicting  the  appearance  of  satellite  photographs,  especially  those 
over  large  areas  such  as  ATS  photographs.  In  the  tropics  and  subtropics 
where  conventional  weather  models  (such  as  the  wave  cyclone;  are  not 
applicable,  the  compensation  model  is  useful. 

For  the  purpose  of  exploinmg  this  model,  the  condition  is  accepted 
that  low  level  flow  under  a  trade  or  other  tropical  inversion  is  very  constant 
and  very  persistent  compared  to  the  upper  level  flow.  Most  of  the  time  the 
inversion  rises  or  falls  slightly  to  allow  compensation  for  upper  air  pressure 
changes  and  thus  leaves  the  lower  flow  undisturbed.  This  process  is  the 
source  of  the  term. 

If  the  inversion  height  over  the  lower  moist  layer  is  near  the 
condensation  level,  or  if  cloud  heights  can  be  determined,  then  the  changes 
in  height  of  the  inversion  will  be  indicated  by  clouds. 

Some  cloud  patterns  to  be  expected  with  traditional  flows  are 
shown  in  Figures  1-4.  Easterlies  with  no  upper  air  disturbance  and 
light  upper  wind  coulci  have  the  inversion  getting  higher  to  the  north  as 
indicated  in  Figure  1.  Arrows  show  upper  wind,  barbs  show  surface 
wind.  The  pattern  with  an  upper  air  trough  would  indicate  that  the 
inversion  rises  slightly  to  keep  the  trough  from  showing  up  in  the  surface 
flow.  This  is  shown  in  F'.gure  2,  which  gives  an  indication  of  the  expected 
cloua  pattern. 
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Figure  1.  Easterlies  with  no  upper  air  disturbance  and  light  upper  wind 
will  have  higher  inversion  in  the  north.  High  inversion  height  is  shown 
by  thicker  cloud 


Figure  l  .  Upper  air  trough  causes  inversion  to  rise  with  little  effect 
on  surface  flow.  High  inversion  height  is  shown  by  thicker  clouds. 


Figure  3.  Upper  air  ridge  causes  Inversion  to  lower  suppressed  clouds. 
Higher  Inversion  height  is  shown  by  thicker  clouds. 
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Figure  4.  Upper  all  low  cel!  causes  Inversion  similar  to  trough,  with  modlfi 
cation  as  shown.  Higher  Inversion  heights  are  shown  by  thicker  cioubs. 


An  upper  aLr  ridge  wouid  cause  the  inversion  to  get  lower  and 


would  suppress  the  clouds  as  in  Figure  3.  An  upper  air  low  cell  would 
be  similar  to  a  trough,  but  with  some  modification,  as  shown  in  Figure  4. 
There  would  be  a  similar  modification  for  a  high  cell.  These  patterns 
would  be  particularly  meaningful  if  the  satellite  photographs  were  watched 
daily  and  changes  noted.  For  example,  Figure  5b  would  indicate  a  trough 
and  Figure  5c  would  indicate  a  ridge  over  the  area  which  originally  looked 
like  5a. 

To  summarize,  compensation  is  a  slight  adjustment  of  the  inversion 
height  to  compensate  for  upper  air  pressure  changes  and  leave  the  flow 
below  the  inversion  relatively  undisturbed.  Response  In  the  form  of  height 
change  is  small,  and  vorticity  changes  in  the  lower  layer  are  usually 
negligible. 
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f  igure  S .  Prediction  ot  inversion  from  satellite  photographs,  b.  would 
indicate  a  trough,  cr  c.  would  indicate  a  ridge,  over  area  previously  like 
a.  .  That  is,  higher  inversion  heights  are  shown  by  thicker  clouds. 


C.  Moist  Qvercompensation 

When  clouds  form  in  a  moist  layer  forming  an  inversion  below 
a  dry  layer  there  is  a  built-in  mechanism  for  weakening  the  inversion. 

As  the  air  rises  the  saturated  lower  air  cools  with  the  moist  adiabatic 
lapse  rate.  If  the  air  continues  to  rise  the  inversion  will  eventually 
disappear.  This  mechanism,  coupled  with  the  compensation  mechanism, 
leads  to  moist  overcompensation.  When  the  height  required  for  compen¬ 
sation  is  so  great  that  the  inversion  must  exceed  the  level  at  which  the 
inversion  is  dissipated,  then  moist  overcompensation  results. 

Formation  of  clouds  means  that  increased  inversion  height 
becomes  less  and  less  effective  in  causing  compensating  pressure 
changes.  Decreasing  pressures  aloft  can  only  be  compensated  by  extra¬ 
ordinarily  high  inversion  heights,  resulting  in  thick  clouds  and  possible 
measurable  changes  in  the  vertical  component  of  the  vorticity.  In 
extreme  cases  all  stability  is  eliminated.  It  finally  becomes  impossible 
for  the  surface  pressure  pattern  to  differ  from  the  upper  air  pressure  pattern 
and  the  surface  flow  comes  under  the  influence  of  strong  pressure 
gradients  from  aloft. 

This  phenomenon  can  best  be  illustrated  by  describing  the  flow 
in  the  easterlies  under  a  deepening  upper  air  trough.  In  Figures  6,  7,  8, 

9,  and  10  we  have  plotted  contours  of  the  inversion  height,  upper  air 
winds  and  surface  winds. 
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Figure  6.  Schematic  represen¬ 
tation  of  undisturbed  subtrop¬ 
ical  easterly  flow  underlying 
upper  level  westerlies.  Upper 
level  westerlies  superimposed 
on  inversion  height  contour 
lines . 


Figure  7.  Illustration  of  com¬ 
pensation  by  upward  adjust¬ 
ment  of  Inversion  height  for  a 
disturbance  in  the  upper 
westerlies  so  that  low  level 
easterlies  are  maintained 


Figure  8  Illustration  of  the 
case  where  compensation 
cannot  maintain  the  low  1. .  1 
easterly  flow,  and  the  hlch  level 
perturbation  Is  reflected  at  the 
surface . 


Figure  9-  Evolution  of  low 
level  flow  ana  cloud  patterns 
subsequent  to  Figure  8.  Inver¬ 
sion  has  been  destroyed  in  clom 
area  In  all  drawings  the  inver¬ 
sion  contours  are  almost  para¬ 
llel  to  the  upper  winds  . 
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rigurr  in.  Large  scale  moist  overcompensation  leads  to  the  upper  air 
pressure  gradient  causing  "geostrophic "  west  winds  at  the  surface. 
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Figure  6  shows  the  straight  westerly  flow  over  the  straight 
easterly  flow,  with  the  inversion  sloping  to  compensate  for  the  dif¬ 
ference  between  fiows.  The  low  pressure  in  the  upper  trough  in 
Figure  7  is  compensated  by  the  increase  in  inversion  height  so  the 
east  winds  continue  in  the  lower  levels  undisturbed.  In  Figure  8 
the  trough  has  become  so  strong  that  the  compensating  effect  has 
lifted  the  air  so  high  the  clouds  formed  have  reduced  gravitational 
stability  to  almost  nothing.  Finally,  it  is  no  longer  possible  to 
compensate  and  the  upper  air  pressure  gradient  cannot  be  shielded 
from  the  surface  air.  Figure  9  Illustrates  a  strong  south  wind  resulting 
from  acceleration  of  the  surface  air  by  the  upper  air  pressure  grad¬ 
ient  that  accompanies  west  winds  aloft. 

If  this  event  occurs  on  a  large  enough  scale  so  that  geo- 
strophic  balance  is  expected,  then  we  have  west  winds  at  the  ground 
as  shown  in  Figure  10  . 

On  a  smaller  space  scale  or  a  shorter  time  scale  it  would  be 
expected  that  the  air  moving  into  the  cloudy  area  under  west  winds 
aloft  would  move  to  the  north  (Figure9).  This  northward  moving  air 
is  not  necessarily  confined  to  the  surface  layers,  so  it  can  continue 
to  rise  and  become  a  plume  of  clouds  at  high  levels  streaking  from 
southwest  to  northeast. 
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D.  Moist  Air  Injection 

Occasionally  the  southerly  {or  subtropical)  core  of  the  jet 
stream  undulates  equatorward  over  the  moist  easterlies,  and  causes 
moist  overcompensation.  When  the  ability  for  compensation  is  exhaus¬ 
ted  the  warm  moist  air  rises  and  pushes  up  the  tropopause,  causing  an 
increased  tropopause  height  over  a  large  area.  This  increased  height 
starts  two  waves  on  the  tropopause,  one  streaming  off  to  the  southeast  and 
one  streaming  off  to  the  northeast.  The  wave  streaming  off  to  the  south¬ 
east  is  lost  in  the  equatorial  cloudiness  but  the  one  moving  off  to  the 
northeast  becomes  a  streak  of  moist  air  and  cloud--and  highly  visible  on 
satellite  photographs  {see  streak  off  west  coast  of  Baja  California  In 
Figure  ZO). 

The  northwest  side  of  this  injection  can  be  in  geostrophic  balance 
and  may  not  contribute  to  the  northward  advance  of  the  wave,  but  the 
southeast  side  cannot  be  in  geostrophic  balance  with  southwest  winds 
so  the  advance  of  the  tropopause  wave  would  be  expected  to  have  a  peak 
behind  it.  (Figure  11). 
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1, .  Theory  of  Moisture  Influx 


The  source  of  moisture  in  the  air  is  nearly  always  considered  to 
be  near  the  ground.  In  order  to  simulate  this  concept  (which  is  also 
backed  by  adequate  observation)  we  assume  that  there  is  a  layer  near  the 
ground  with  gravitational  stability  and  a  high  relative  humidity  (50-70$). 
When  this  moist  layer  is  lifted  enough  that  the  relative  humidity  becomes 
100$,  clouds  form,  and  the  lifting  lapse  rate  in  the  low  level  air  becomes 
the  saturated  adiabatic.  Then  if  the  upper  air  is  lifted,  following  the  dry 
adiabatic  rate,  the  net  result  is  a  decrease  in  gravitational  stability. 
Repetition  of  this  process  leads  to  less  and  less  stability,  enabling  the 
clouds  to  move  very  rapidly  up  to  the  tropopause  where  the  stability 
again  becomes  very  large. 

The  particular  formulation  of  the  problem  developed  here  is  based 
on  gravitational  stability  in  the  rising  air.  This  was  done  primarily  for 
convenience  of  computing,  but  there  is  some  basis  for  it  in  certain  large 
circulation  patterns  where  complete  instability  (warm  moist  air  under  cool 
ary  air)  is  not  likely.  A  large  ,  absolutely  flat,  unstable  region  would  allow 
cooler  dry  air  to  fall  into  the  moist  air,  and  the  resulting  cooler  moist  air 
might  then  be  stable.  See  Figure  12.  The  arrangement  in  (a)  is  not  likely. 
But  if  the  large  flat  area  of  cold  air  is  unstable,  so  that  the  cold  air  intrudes 
downward,  it  will  cause  a  picture  similar  to  (b),  a  stage  (perhaps  not  ob¬ 
served)  in  the  mixing  of  cool  upper  air  into  a  large  area  of  saturated  moist 
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air.  The  actual  picture  Ls  probably  more  like  (c),  which  implies  mixing  of 
the  rising  warm  air  and  dry  air  above,  with  the  cool  mixed  air  flowing 


downward  on  the  edges. 


F.  Introduction  of  Moisture  Into  the  Prediction  Model 

The  study  of  the  effect  of  moisture  on  the  flow  is  carried  out  with 

the  following  system  of  equations.  We  define  the  symbols  below; 

9  —  potential  temperature  below  inversion  in  degrees  Celcius 

9'  —  potential  temperature  above  inversion  in  degrees  Celcius 
A 9  —  difference  In  temperature  across  the  Inversion  at  conden¬ 
sation  level  In  Celcius  degrees 

hc  —  the  condensation  level  in  miles 

H  —  height  of  the  tropopause  (approx.  40,000  ft,  )in  miles 
F i  —  upper  constant  of  stability  in  mi/hr2 
h^  —  maximum  gravitational  stability  height  in  miles 
h  —  height  of  inversion  above  sea  level  in  miles 
hg  —  height  of  ground  surface  above  sea  level  in  miles 
Compute:  hm 

Note  h  here  is  the  maximum  height  for  gravitational  stability;  the 
m 

low  level  air  is  cooler  than  the  upper  level  air  when  h  is  less  than  h 

nr 

From  A6  -  7.2  (h  -  h  )  =  0  obtain  h_  = AS/7. 2  +h  .  The  input  para- 
m  c  m  c 

meters  are  assumed  to  satisfy  the  conditions; 

AS  2  0  and  H  2  h  . 

m 

Figure  13  Illustrates  temperature  distribution.  The  conversion 
factor  used  Ls  1.8.  There  are  four  cases; 

A0 

Case  I  h  <  h  V  =  -r-g 

c  G 
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In  these  equations  u,v  are  wLnd  components  in  the  lower  air. 

u'  ,v'  are  the  wind  components  in  the  upper  air.  h^  is  the  height  of  the 

ground  surface  above  sea  level;  k  is  the  drag  coefficient  of  the  ground 

Cloud  formation  is  associated  with  a  higher  and  a  weaker  inversion. 

This  is  illustrated  by  Case  II,  ( h  £  h  <  h  )  with  Equation  (1)  . 

cm 

When  the  inversion  height  reaches  h  ,  assume  that  gamma  is 

m 

constant  in  the  region  between  the  maximum  height  for  gravitational 
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stability,  h  ,  and  the  tropopause  height  H,  and  we  have  Case  III, 
m 

(h  £  h  ^  H)  with  Equation  (1).  This  small  amount  of  stability  takes 
m 

into  account  the  cooling  and  mixing  of  the  moist  air  with  the  upper  air 
by  maintaining  the  stability. 

When  the  surface  air  reaches  the  tropopause  we  have  a  very 
stable  deep  layer  that  follows  the  equations  In  Case  IV,  (h  >H)  with 
Equation  (1)  but  now  with  u'  and  v'  =  0. 

In  each  of  the  above  cases  the  inversion  level  winds  u'  and 
v'  are  assumed  to  be  functions  of  the  inversion  height  h  .  The 
functions  u'(h)  and  v'(h)  are  defined  so  as  to  approximate  the  ambient 
thermal  winds  between  the  low  level  inversion  and  the  tropopause.  The 
functions  u'  and  v'  could  be  defined  as  u'(x,  y,  t,  h)  and 
v'(x,  y,  t,  h)..  For  this  particular  program  they  were  defined  as  u'(h) 
ana  v'(h). 
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II.  THE  SYNOPTIC  SITUATION  FOR  THE  PREDICTION  EXPERIMENT 


Keith  Veigas  ana  Pieter  Feteris 

A.  Procedure 

The  map  and  ATS  photo  series  for  6-10  April  1967  were  selected 
for  initial  prediction  experiments  employing  the  model  described  in  sec¬ 
tion  I  of  this  report.  As  all  meteorological  maps  are  marked  Z  for  GMT, 
this  notation  has  been  used  for  cons istency  and  brevity. 

Figures  14,  15,  and  16  for  7  April  1967  illustrate  the  initial  condi¬ 
tions  for  the  first  prediction  experiment.  Figure  14  shows  the  surface 
chart  for  0000Z,  and  the  associated  ATS  photocommencing  approximately  2 
hours  prior  at  around  22  17Z  6  April. 

The  shaded  area  on  the  surface  chart  (Figure  14)  is  the  prediction 
area.  The  low  level  flow  is  predominately  easterly  and  the  ATS  photo 
coupled  with  available  surface  observations  indicate  only  very  widely 
scattered  areas  of  convective  activity. 

Circulation  in  the  middle  and  upper  troposphere  at  CC00Z  7  April 
are  shown  on  Figures  15  and  16.  There  is  an  upper  level  vortex  centered 
over  the  western  boundary  of  the  prediction  area,  with  southwesterly  winds 
at  500  mb  over  the  low  level  easterlies  and  predominacely  westerly  flow 
in  the  upper  troposphere. 

Approximately  24  hours  later  the  ATS  photo  commencing  at  2  340Z  7 
April  and  reproduced  in  Figure  17  indicates  a  considerable  difference 
in  the  cloudiness  and  associated  moisture  distributions.  Also  included 
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Figure  17.  Surface  chart  for  OOOOZ  8April  and  ATS  photo  for  2  340Z  7  April 
1967. 
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in  Figure  17  is  the  surface  chart  for  OOOOZ  8  April.  A  plot  of  conven¬ 
tional  cloud  observations  on  this  day  indicates  that  the  cloudy  area 
labeled  A  on  the  ATS  photo  is  marked  by  much  more  widespread  deep  con¬ 
vection  than  on  the  previous  day.  The  flow  In  the  lower  layers  to  the 
south  and  east  of  Hawaii  Is  still  predominantly  from  the  east.  In  the 
mid-troposphere  it  can  be  noted  that  a  cyclonic  vortex  remains  In  the 
western  section  of  the  region  of  Interest,  and  the  associated  southerly 
flow  that  has  developed  In  the  mld-tropospherlc  regions  (Figures  18 
and  19)  is  still  not  reflected  In  the  lower  layers  of  the  atmosphere. 
Qualitatively,  according  to  the  arguments  of  Freeman  (Section  I),  this 
should  lead  to  an  increase  in  the  depth  of  the  zone  of  the  easterlies  from 
Hawaii  toward  the  center  of  the  mid-tropospheric  vortex.  An  Increase  In  the 
depth  of  the  easterlies  implies  upward  vertical  motion  and  eventual  release 
of  convective  instability.  This  can  lead  to  condensation  In  the  moist  lower 
layer , leading  to  the  development  of  vigorous  convective  action  which 
will  "break"  the  cap  of  the  inversion  and  the  development  of  moist  areas 
through  rather  deep  layers  of  the  troposphere.  A  deep  moist  layer  can 
subsequently  be  Injected  into  midlatitude  weather  systems.  An  evolving 
pattern  of  this  type  is  indicated  in  Figures  20  and  21  of  the  ATS  photos 
late  on  the  8th  and  9th  of  April. 
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Figure  18.  500  mb  stramline  analysis  chart  for  0000Z  8  April  1967. 
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Figure  19  2  50  rr.b  streaml 


III.  RESULTS  OF  THE  NUMERICAL  PREDICTION  EXPERIMENTS 


Keith  Veigas  and  Sam  Rosenberg 

As  Indicated  In  section  II  of  this  report  the  weather  satellite 
sequence  commencing  OOOOZ  7  April  1967  was  selected  for  our  initial 
prediction  experiments. 

The  important  features  of  the  synoptic  situation,  as  revealed  by 
the  ATS  photograph  and  associated  synoptic  charts  shown  In  Figures 
14,  15  and  16  ,  are  a  broad  belt  of  low  level  easterly  winds,  only  widely 
scattered  convective  activity,  and  a  middle  and  upper  tropospheric  dis¬ 
turbance  at  the  western  border  of  the  prediction  area. 

Two  prediction  experiments  were  carried  out.  The  height  of  the 
zone  of  the  easterlies  was  assumed  to  be  everywhere  equal  to  one  mile; 
In  the  first  experiment  £8  ,  strength  of  the  inversion,  was  designated 
3°C  and  in  the  second  6°C. 

The  upper  boundary  conditions  and  those  at  the  periphery  of  the 
grid  were  held  constant,  which  limited  the  prediction  period  to  six 
hours . 

Estimates  of  the  T0  winds  were  based  on  analyses  prepared 
by  the  ESSA  Weather  Bureau  In  Honolulu.  Analyses  were  available  for 
the  surface,  700,  500,  and  250  mb  levels.  The  components  of  the  T0 
low  level  wind  field  are  shown  in  Figures  22  and  23. 
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The  components  of  the  low  level  wind  field  at  t0  +  6  hours  are  shown 
in  Figures  24  and  25.  The  u  component  exhibited  only  small  changes, 
but  a  rather  dramatic  change  in  the  v  component  from  nearly  zero  to  a 
strong  southerly  component  was  predicted.  (This  predicted  change  was 
realized  (see  Figure  17);  however,  the  zone  of  Increase  of  southerly 


winds  was  northwest  of  the  predicted  area  of  Increase.) 

The  predicted  change  in  height  of  the  layer  of  the  easterlies  Is  shown 
in  Figure  26  and  there  is  a  good  correspondence  between  the  area  of 
implied  upward  vertical  motion  and  the  developing  cloud  areas  as  depicted 
by  a  subsequent  ATS  photograph  (Figure  17)  for  2340Z  7  April  1967.  For 
example,  see  cloud  areas  labeled  6,0,0,  on  the  prognostic  h  chart 
and  the  ATS  photograph.  In  each  case  the  upward  motion  implied  by  the 
h  prognostic  chart  corresponds  quite  well  with  the  developing  cloud  areas. 

Experiment  II  with  6 6  set  at  6°C  produced  results  very  stmilar  to 
experiment  I  illustrated  in  Figures  24,  25,  and  26.  The  predicted  southerly 
components  of  the  low  level  wind  and  the  increases  in  the  height  of  the 
low  level  patterns  were  slightly  higher,  but  It  can  be  concluded  that 
the  model  is  not  highly  sensitive  to  the  specification  of  66.  This  Is  an 
important  result  in  that  only  gross  estimates  of  60  can  be  expected  In 
the  data  sparse  regions. 

It  is  felt  that  future  diagnostic-prediction  studies,  which  are  beyond 
the  scope  of  this  study,  should  include: 

1 .  Introduction  of  time  dependent  boundary  conditions  from  mld- 
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Figure  26.  Map  of  predicted  height  (nautical  miles)  of  the  layer  of  the 
easterlies  0600Z  7  April  1967. 


latitude  numerical  predictions  where  possible. 

Introduction  of  effects  of  the  release  of  latent  heat  of 
condensation  from  convective  scale  processes  on  the 
synoptic  scale  generation  of  kinetic  energy. 

Introduction  of  the  flux  of  generated  kinetic  energy  from 
low  latitudes  into  mio-latitude  prediction  models. 

Designing  data  archiving  programs  for  satellite  and  con¬ 
ventional  data ,  so  that  statistical  "tuning"  of  this  and  other 
theoretical  models  is  possible. 

Utilizing  special  data  collections,  such  as  those 
anticipated  in  Project  Bomex,  to  carry  out  more  definitive 
prediction-diagnostic  experiments. 


IV,  ON  THE  FUTURE  USE  OF  ATS  TYPE  DATA 


IN  NUMERICAL  PREDICTION  MODELS 
Keith  Veigas 


A.  Case  Study 

In  the  course  of  examining  a  series  of  weather  maps  and  associ¬ 
ated  ATS  photos  in  order  to  select  test  cases  for  the  application  of  a 
dynamical  prediction  model  for  this  study,  several  interesting  situations 
were  encountered  which  supply  clues  to  ways  in  which  ATS  type  information 
may  be  introduced  into  present  and  future  "wet"  numerical  prediction  models. 
In  addition  to  the  map  series  described  in  Section  I,  the  series  extending 
from  2 6  to  29  March  1967  is  notewortr 

Figure  27  shows  the  surface  analysis  for  0000Z  23  March  19o7  and 
the  associated  ATS  photo  for  2252Z  27  March  1967.  Note  the  cloud  area 
B  on  the  satellite  photo  and  the  surface  chart.  It  is  obvious  from  the  con¬ 
ventional  cloud  observations  that  the  cloud  area  represents  a  deep  layer 
of  moisture  that  has  developed  in  the  warm  sector  of  the  wave  cyclone 
centered  at  about  39°N,  176°W  longitude.  The  cyclone  had  been  moving 
in  a  rather  normal  and  consistent  track,  as  indicated  by  the  six  hourly 
positions  shown  in  Figure  27. 

Figure  28  illustrates  some  racher  dramatic  changes  in  the  cyclonic 
system.  The  surface  analysis  in  Figure  28  is  for  1800Z  28  March  and  the 
ATS  photo  is  for  2137Z  28  March  1967.  Judging  from  a  careful  analysis  of 
the  sequence  of  hourly  surface  charts  there  is  good  observational  evidence 
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Figure  27.  Surface  chart  for  OOOOZ  28  March  and  ATS  photo  for  22  52Z 
29  March  1967. 


that,  the  low  cited  in  the  previous  figure  split  into  two  centers  with  the 
northernmost  center  moving  to  the  ENE  where  it  dissipated  before  the 
1800Z  28  March  map  time.  The  southern  center  moved  SE  and  deepened. 

An  additional  cyclone  formed  under  the  center  of  the  deep  moist  layer  indi¬ 
cated  in  Figure  27.  Upon  this  new  development  the  axis  of  these  cy¬ 
clones  rotated  in  a  counterclockwise  manner. 

In  studying  the  sequence  of  events  as  depicted  by  the  ATS  and  other 
satellite  photos  quite  a  different  picture  evolves.  There  is  a  noticeable 
shift  to  the  south  of  the  cloud  mass  associated  with  the  cyclone  centered 
about  40°N  latitude  at  OCOOZ  28  March,  and  the  cloud  area  denoted  by  B 
in  Figure  27  was  injected  into  the  midlatitude  weather  system.  Thus  it 
seems  that  utilizing  satellite  photos  would  likely  have  led  to  a  better 
surface  prediction  in  that  it  did  indicate  the  tendency  for  the  midlatitude 
weather  system  to  shift  to  a  more  southerly  course,  and  perhaps  more  Im¬ 
portantly,  the  forecast  of  weather  e.'cn  ?nts  could  likely  have  been  Improved 
with  the  knowledge  of  tne  spatial  extent  and  direction  of  motion  of  a  deep 
moist  layer  advancing  into  the  weather  system  from  the  south. 

B  Use  of  Satellite  Infer. nation  in  Numerical  Weather  Prediction 

While  useful  information,  i.e.  circulation  parameters,  maybe 
inferred  from  ATS  ty;  .?  photos,  additional  useful  information  may  be  In¬ 
tro  limed  in  NWP  cloud  and  precipitation  models  by  the  use  of  future  ATS 
satellites.  There  is  also  increasing  evidence  that  convective  scale 
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release  of  latent  heat  In  the  tropics  and  subtropics  must,  in  some  manner, 
be  taken  into  account  in  future  prediction  models. 

The  recent  studies  of  hurricane  development  by  Kuo  (1965)  and 
Miller  (1969)  clearly  show  that  predicted  behavior  of  forecast  models  de¬ 
pend  very  critically  on  the  moisture  content  of  the  air.  Blending  of  a 
variety  of  types  of  available  information  will  be  required  to  obtain  the 
best  possible  estimates  of  moisture  distributions  over  data  sparse  regions 
of  the  globe  for  input  to  future  operation  tropical  and  subtropical  pre¬ 
diction  models.  It  seems  obvious  from  the  ATS  photos  studied  in  the 
course  of  preparing  for  the  numerical  prediction  experiments  outlined  in 
sections  I  and  III  of  this  report  that  ATS  information  can  be  extremely 
useful  in  delineating  the  spatial  extent  of  vertical  moisture  distributions 
detected  at  points  from  conventional  observation  sources. 

Prediction  experiments  should  be  designed  to  seek  means  of 
effectively  incorporating  this  new  source  of  Information  into  operational  and 
research  prediction  models,  and  further,  to  Isolate  the  impact  of  improved 
initial  moisture  specification  on  the  prediction  of  clouds,  precipitation  and 
circulation. 

In  addition  it  seems  timely  to  suggest  that  the  development  of  methods 
to  combine  various  sources  of  humidity  Information  should  be  accelerated. 
Chisholm,  Ball,  Velgas,  and  Luty,  (1968)  and  Ball  and  Velgas,  (1968)  have 
shown  that  it  is  feasible  to  combine  several  types  of  humidity  information 
to  estimate  free  atmosphere  humidity.  They  have  further  demonstrated  that 
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these  estimates  lead  to  improved  objective  analysis  of  humidity  dis¬ 
tributions  in  data  sparse  regions. 

The  coupling  of  conventional  and  satellite  information  is  a  difficult 
problem,  and  much  care  must  be  exercised  in  developing  techniques  for 
inclusion  of  this  new  Information  into  forecasting  systems.  This  new 
information  should  be  an  invaluable  adjunct  to  presently  available  infor¬ 
mation,  as  illustrated  by  Figure  29.  The  USAF  AWS  Global  Weather  Center 
kindly  furnished  us  with  their  objective  cloud  analysis  for  the  0000Z  28  March 
1967  study  discussed  above.  Toe  hatched  area  in  Figure  2  9  is  that  area 
where  the  cloud  amount  was  classified  as  indeterminate  by  the  GW C 
analysis  program  because  of  lack  of  data.  With  the  operational  use  of 
ATS  satellites  it  is  hoped  that  these  "silent  areas"  will  cease  to  exist. 

In  summary,  these  case  studies  have  indicated  that; 

1.  The  use  of  satellite  data  to  augment  conventional  observations 
will  lead  to  more  accurate  estimates  of  moisture  parameter  for 
input  to  NWP  models. 

2.  It  became  obvious  tnat  combining  conventional  and  satellite 
data  is  not  a  trivial  problem,  because  the  different  types  of 
information  vary  in  characteristics  and  accuracy. 

3.  Moisture  inputs  plus  satellite-derived  sounding  inputs  should 
yield  much  better  specification  of  "wet"  NWP  model  initial 
conditions . 
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Figure  29.  Hatched  areas  denote  zones  where  conventional  observations 
are  inadequate  to  specify  cloud  cover  on  0000Z  28  March  1  967  . 
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Definitive  prediction  experiments  should  now  be  designed  to 
see  if: 

a.  Improved  cloud  and  precipitation  predictions  can  be  made. 

b.  Existing  operational  end  research  models  can  properly 

accommodate  conversion  of  latent  heat  energy  to  kinetic 
energy  in  our  circulation  predictions. 


REFERENCES 


Armijo,  Larry,  1966:  Statistical  properties  of  radar  echo  patterns  and 
the  radar  echo  process.  Technical  Note  48-NSSL-27.  National 
Engineering  Science  Company,  Houston. 

Ball,  J.  T.  and  K.  W.  Veigas,  1968:  The  analysis  of  upper-level  hu¬ 
midity.  J.  Appl.  Meteorol.  ,  7,  4,  620-625. 

Bomex  Bulletin  No.  4,  prepared  by  The  Bomex  Project  Office,  May  1969. 

Chisholm,  Donald  A.;  John  T.  Ball;  Keith  W.  Veigas,  and  Paul  V.  Lutz, 
1968:  The  diagnosis  of  upper-level  humidity.  J.  Appl.  Meteorol 
7,  4,  613-619. 

Kessler,  Edwin  III;  and  John  A.  Russo,  Jr.  ,  1963;  A  program  for  the 
assembly  and  display  of  radar-echo  distributions,  J.  Appl. 
Meteorol.  .  2  .  5,  October,  582  -593. 

Kuo,  H.  L.  ,  1965;  On  formation  and  intensification  of  tropical  cyclones 
through  latent  heat  release  by  cumulus  convection.  J.  Atmos. 
Sci.  ,11,  1,  40-63. 

Leese,  John  A.  ,  and  Edward  Epstein,  1963;  Application  of  two-dimen¬ 
sional  spectral  analysis  to  the  quantification  of  satellite  cloud 
photographs.  J.  Appl.  Metorol.  ,  2  ,  5,  October,  629-644. 

Miller,  B.  I.  ,  1969:  Experiment  in  forecasting  hurricane  development 
with  real  data.  Nat1 1,  Hurr .  Res .  Proj .  ,  Rept.  No.  8  5. 


47 


■t  m**  itt:  fi* W  '[r  • 


j 


i 

i 

APPENDIX  A 

A  METHOD  OF  FOLLOWING  PATTERN  CLOUDS  ON 
SEQUENTIAL  SATELLITE  PHOTOGRAPHS 
John  C.  Freeman  and  Troxel  Ballou 

I  A.  Discussion 

t 

i 

k 

f  The  recent  spectacular  success  of  the  ATS-1  and  AT3-3  synchro- 

I 

|  nous  satellites  has  made  it  obvious  that  clouds  can  be  photographed 

i 

6 

j  over  large  areas.  This  appendix  outlines  one  numerical  method  for  deter- 

V 

T 

i  mining  motions  of  clouds  from  these  photographs. 

t 

The  current  most  accepted  method  of  following  cloud  motions  for 
|  the  purpose  of  making  research  studies  is  to  make  a  movie  from  the  ATS 

photographs  and  follow  the  cloud  motion  in  the  movie.  The  most  useful 
of  these  movies  have  been  made  by  Fujita.  Similar  devices  involving 
flickers  of  these  photographs  have  been  developed  in  the  National  Environ- 

} 

» 

mental  Satellite  Center. 

!  The  problem  of  determining  the  motion  of  clouds  is  part  of  the 

satellite  photograph  modeling  problem.  The  method  described  here  has  not 
;  been  investigated  sufficiently  to  give  any  idea  of  its  feasibility  or  its  ex- 

I  pense  but  it  seems  worthy  of  including  in  this  appendix  as  a  guide  to  future 

i 

discussion.  The  method  requires  the  following  equipment: 
a.  A  light  intensity  intergrator  with  digital  output 
j  b.  An  electronic  computer 

i  i 
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c.  Image  reproducing  equipment 

If  a  forecast  system  is  being  developed  the  measured  cloud 
speeds  would  be  used  to  move  the  cloud  patterns.  This  could  be  done 
in  a  selective  way  by  printing  the  matrices  for  the  speed  maxima.  An 
example  of  how  this  works  is  shown. 

This  method  of  pattern  following  and  especially  recognition  of 
patterns  has  been  advanced  in  closely  related  forms  by  Kessler  and  Russo 
{  1  96 3 j ,  Leese  and  Epstein  (1963)  and  Armijo  (1966).  However  the  method 
described  here  is  simpler  and  has  a  more  topical  objective  than  any  of 
those. 

B.  Procedure  for  Analysis  of  Cloud  Movement  on  Satellite  Radar 
Cloud  Maps 

1.  Select  two  maps  of  the  same  area  at  times  tj  and  t£. 

2.  Select  a  large  grid  mesh  (e.g.  10-29  degrees). 

3.  Divide  it  into  gridpomts  no  more  than  1/10  the  size  of  the 
large  grid. 

4.  Digitize  the  data  at  each  gridpoint  of  both  maps  into  binary 
digits  by  some  criterion  (such  as  "clouds/no  clouds"  present 
at  that  gridpoint). 

5.  Select  a  small  grid  mesh  for  correlating  (see  fa).  A  reasonable 
size  would  be  five  times  the  gridpoint  spacing;  an  odd  number 


greater  than  three  is  recommended. 


6.  For  each  gridpoint  away  from,  the  boundary  (for  points 

near  the  boundary,  the  overlay  extends  past  the  boundary) 
correlate  the  data  at  surrounding  points  in  this  manner: 


let  k  range  from  i  -  1  to  i  +  1 


let  i  range  from  j  -  1  to  )  +  1 

for  all  combinations  of  t^  overlaid  by  tj.., 

define  the  number  s^  as  the  sum  of  the  product 

ti  t?  for  all  m,n  under  the  sample  grid, 
mn  mn 

7.  Of  these  nine  sums,  the  largest  gives  the  magnitude  of  a 
vector  with  direction  numbers  (k  -  i,  l  -  j)  which  in¬ 


dicates  the  trend  of  cloud  motion  at  ti.  . 

U 


Example; 


Figure  30a  shows  a  15°  square  grid  with  qridpoints  1°  apart. 

Blank  squares  contain  zeroes,  other  squares,  ones,  x  indicates  the 
center.  Select  a  sampling  grid  5°  square  (Figure  30b).  At  the  point 
(t,j)  overlay  tj  onto  t^  (Figure  30c).  For  the  35  squares  that  over¬ 
lap,  take  the  product  of  each  top  square  with  the  one  directly  below  it 
(Figure  31).  Sum  these  35  numbers  and  call  the  sum  .  In  this 
example,  s^  is  4.  Repeat  this  procedure  eight  more  times,  super¬ 
imposing  on  each  of  the  eight  points  surrounding  t,  .  Note 
ij  kZ 

that  to  do  this  it  is  necessary  to  have  values  for  an  extra  gridpoint  all 
the  way  around  the  section  of  tz  grid  which  is  overlaid. 
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(C) 


Figure  30  Binary  digitized  map  of  satellit'-  radar  cloud  photographs, 
(a)  Full  grid,  with  binary  digitized  data,  (b)  Sampling  grid  from  t! 


centered  at  t, 
centered  at  tt 


(c)  Urid  around  t?.  overlaid  by  sampling  square 
k* 


Figure  3  1  Section  from  Figure  30(c).  Since  there  are  four  positions  con¬ 
taining  ones  at  coincident  gridpoints  above  and  below,  the  sum  a  ,  is  4. 
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Define  a  vector  which  points  from  t4  to  the  maximum  or 

center  of  gravity  of  these  nine  sums,  and  has  a  magnitude  equal  to 

that  maximum:  e.g.  ,  if  the  maximum  sum  was  at  t2.  .  ,  the  direction 

1  +  i  /  } “  1 

numbers  of  the  vector  would  be  (1,-1),  and  its  magnitude  s.  .  .  . 

I  T"  1  ,  J  **  1 

This  vector  is  the  primary  gross-scale  estimator  of  cloud  movement. 

Figures  32-34  show  an  actual  example.  The  data  were  from  two 

ESSA  digitized  mosaic  maps  for  29  and  30  May,  1969  at  2000Z  (Figures  35 

and  36)  .  The  original  grid  was  10°  square,  the  gridpoint  spacing  1°,  and 

the  sampling  grid  size  4°  (5  would  have  been  too  large  for  this  small 

a  grid;  4  was  small  enough,  but  difficult  to  manipulate).  Figure  32 

shows  the  grid  after  all  the  sums  have  been  taken.  The  nine  numbers 

in  each  square  (blanks  indicate  zeroes)  are  the  magnitude  of  the  vector 

pointing  in  that  direction  from  the  center  of  the  square.  Fcr  instance, 

in  the  square  at  7  5°W  24°N,  the  5  means  that  when  ti_.  , ,  was  over- 

75,24 

laid  onto  t2  ,  there  were  5  squares  that  matched  (with  ones).  With 

i  D  t  CD 

*!75  24  over^ai°  on  *^74  25  '  t*lere  were  4  matches. 

In  Figure  3  3  ,  the  highest  number  and  the  primary  trend  has  been 
removed  to  show  secondary  trends. 

In  Figure  34,  the  highest  number  and  the  secondary  trend  in 
Figure  33  has  been  removed  to  show  tertiary  trends. 

The  variance  in  the  direction  of  the  vectors  is  explained  by  the  fact 
that  each  is  weighted  according  to  the  distance  and  magnitude  of  the 


numbers  in  that  square. 
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APPENDIX  B 

SATELLITE  OBSERVED  CLOUDS  AS  STREAK  LINES  AND  THE 
USE  OF  STREAK  LINES  AS  A  COORDINATE  SYSTEM 
John  C.  Freeman 

ABSTRACT 

The  primary  atmospheric  model  emphasized  by  the  satellite 
Is  the  streak  line  of  high  clouds,  that  originate  in  equatorial  regions 
and  carries  a'  jet  stream  of  cloudy  air  all  the  way  up  into  temperate 
latitudes.  This  streak  line  carries  a  companion  dry  streak  -  (to  the 
north)  that  indicates  that  a  systematic  lineal  disturbance  of  the  atmos¬ 
phere  is  closely  related  to  these  streak  lines. 

The  carrying  of  very  narrow  zones  over  great  distances  and 
through  latitude  changes  indicates  that  there  is  some  scientific  infor¬ 
mation  to  be  gained  by  setting  up  a  coordinate  system  based  on  the 
streak  lines.  Such  a  coordinate  system  could  have  the  advantage  of 
being  visible  on  a  satellite  photograph  and  would  make  the  following 
of  physical  phenomena  by  means  of  satellite  data  much  easier  than  it 
is  now. 

C 

The  purpose  of  this  appendix  is  to  develop  a  workable  system  of 

following  quantitative  Information  in  a  streak  line.  The  manuscript  is  avail¬ 
able  from  the  Institute  for  Storm  Research. 

-■This  lias  been  discovered  by  NIMBUS  infrared  data,  found  and  com¬ 
municated  to  us  by  V.  J.  Oliver  (personal  communication,  1%9). 
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13  A  DSTN  AC  T 


Two  models  of  atmospheric  flow,  "compensation"  and  "moist  over-compensa¬ 
tion"  are  found  to  be  useful  in  the  study  of  tropical  weather  in  conventional  and 
satellite  observations.  A  numerical  weather  prediction  model  simulating  important 
physical  processes  in  the  low  latitudes  is  developed  to  test  the  feasibility  of  using 
ATS  information  to  infer  future  weather  distribution.  Diagnostic  studies  using  ATS 
films  and  conventional  synoptic  weather  maps  describe  the  weather. situations  in  the 
Central  North  Pacific  in  terms  of  the  theoretical  models.  Results  of  the  prediction 
experiment  and  future  use  of  ATS  type  data  in  numerical  prediction  models  are  dis¬ 
cussed  and  application  recommended.  Appendix  A  gives  a  description  of  one  method 
of  following  cloud  motion  in  satellite  photographs.  ' 
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